Low- Temperature Heat Transport in the Low-Dimensional Quantum Magnet 

NiCl2-4SC(NH2)2 
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We report a study of the low-temperature thermal conductivity of NiCl2-4SC(NH2)2, which is a 
spin-1 chain system exhibiting the magnon Bose-Einstein condensation (BEC) in magnetic field. It 
is found that the low-T thermal conductivity along the spin-chain direction shows strong anomalies 
at the lower and upper critical fields of the magnon BEC state. In this state, magnons act mainly 
as phonon scatterers at relatively high temperature, but change their role to heat carriers upon 
temperature approaching zero. The result demonstrates a large thermal conductivity in the magnon 
BEC state and points to a direct analog between the magnon BEC and the conventional one. 

PACS numbers: 66.70.-f, 75.47.-m, 75.50.-y 



Bose-Einstein condensation (BEC) denotes the forma- 
tion of a collective quantum ground state of identical par- 
ticles obeying Bose statistics. This fascinating state of 
matter is well established for liquid ''He and ultra-cold 
alkali atoms. It turns out that a form of BEC can also 
be observed in quantum magnets P, 0, H, 0, 01 1 i'^ which 
the density of magnons (bosons) can be tuned by an ex- 
ternal magnetic field (playing the role of chemical poten- 
tial) . Recently, this so-called magnon BEC state has been 
experimentally realized in a growing number of dimer- 
ized spin-1/2 systems, such as the three-dimensional 
system TlCuCla ^,7,], the quasi- two-dimensional sys- 
tem BaCuSi206 [8|, 1^, and the spin- ladder compounds 
(CH3)2CHNH3CuCl3 and (C5Hi2N)2CuBr4 [13, [uj etc. 
A common characteristic of these quantum magnets is 
the spin-gapped ground state at zero field. The external 
magnetic field can close the spin gap and lead to a long- 
range antiferromagnetic (AF) order when the Zeeman en- 
ergy overcomes the gap between the singlet ground state 
and the excited triplet states. It is notable that some 
alternative theories have also been proposed to explain 
the field-induced AF ordered state in quantum magnets 
and challenges the the validity of magnon BEC scenario 
for these systems 12, 13]. To firmly establish the BEC of 
magnons and to develop a deeper understanding of this 
novel state of matter, it would be desirable to look for 
obvious analogs in the basic physical properties between 
the magnon BEC and the conventional one. One of the 
outstanding properties of the superfluid ''He is the ex- 
tremely large thermal conductivity (k) 3, which is 
well understood using the two-fluid model. It is natural 
to ask whether the thermal conductivity (by magnons) 
in the magnon BEC state behaves similarly as that in 
the superfliud ''He. The first experimental exploration 
done by Kudo et al. [161] did reveal an enhancement 
of thermal conductivity at the magnon BEC transition 
of TlCuCla. However, without carefully studying the 
anisotropic heat transport, it is not clear whether the 
enhancement is caused by the weakening of phonon scat- 
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FIG. 1: (color online) Temperature-field phase diagram of 
DTN {H \\ c) obtained from magnetization, specific heat and 
magnetocaloric effect measurements 17, 18, 19, 20]. PM and 
SP represent the low-field quantum paramagnetic state and 
high-field spin-polarized state, respectively. Magnon BEC is 
a magnetic-field induced AF ordered state, with the lower 
and upper critical fields Hci (~ 2 T) and H^2 (~ 12 T). The 
maximum of field-dependent critical temperature is about 1.2 
K, where Hci and Hc2 merge. Inset: Unit cell of tetragonal 
structure showing Ni and CI atoms. The other atoms are 
omitted for clarity. 



tering or the appearing of magnetic heat carriers. 

The organic compound NiCl2-4SC(NH2)2 [dichloro- 
tetrakis thiourea- nickel (II), abbreviated as DTN] is 
the only quantum spin-1 system, rather than spin-1/2 
dimers, to exhibit the BEC of spin degrees of freedom 
[13, [isl, E, [13, EH, El]. It has a tetragonal crystal 
structure (space group 14) [2^, which satisfies the ax- 
ial spin symmetry requirement for a BEC. The Ni spins 
are strongly coupled along the tetragonal c axis (Fig. 
1), making DTN a system of weakly interacting spin-1 
chains with single-ion anisotropy larger than the intra- 
chain exchange coupling. The anisotropy, intra-chain and 
inter-chain exchange parameters of Ni spins were deter- 
mined to be D = 8.9 K, = 2.2 K and Jab = 0.18 



2 



K [13, [il [H, 01, respectively. It was found that the 
Ni spin triplet is split into a 5*2 = ground state and 
5"^ = ±1 excited states with an anisotropy gap of D, 
which precludes any magnetic order at zero field. When 
a magnetic field is applied along the c axis, the Zeeman 
effect lowers the 5^ = 1 level until it becomes degenerate 
with Sz = ground state at Hd , which is essentially the 
same as that occurs in the spin-1/2 dimerized systems. 
Between Hd and Hc2 and below the maximum Tc{H) ~ 
1.2 K (Fig. 1), the magnetic field induces an AF order or 
a magnon BEC. For magnetic field perpendicular to the 
c axis, however, the Sz = ground state mixes with a lin- 
ear combination of the 5^ = ±1 excited states and there 
is no level crossing with increasing field and therefore no 
magnetic order [l7|. DTN was found to be an ideal sys- 
tem for studying the magnon BEC in the sense that its 
lower and upper critical fields are not high, about Hd ^ 
2 T and Hc2 ^ 12 T, which are easily achievable by the 
common laboratory magnets. In this Letter, we show 
that the low-temperature and high-magnetic-field ther- 
mal conductivity of DTN single crystals indeed demon- 
strates a large thermal conductivity in the magnon BEC 
state. 

High-quality NiCl2-4SC(NIl2)2 single crystals are 
grown from aqueous solution of thiourea and nickel chlo- 
ride [is'l. The typical size of single crystals is (0.5-2) 
X (0.5-2) X (3-4) mm"^. X-ray diffraction indicates that 
the parallelepiped crystals are grown along the c axis (the 
maximum dimension) while the four side surfaces are the 
(110) crystallographic plane. So it is easy to prepare 
samples for the thermal conductivity measurements ei- 
ther along (kc) or perpendicular to (nab) the direction of 
spin chains. The thermal conductivity is measured using 
a conventional steady-state technique and two different 
processes: (i) using a "one heater, two thermometers" 
technique in a ^He refrigerator and a 14 T magnet for 
taking data at temperature regime of 0.3-8 K; (ii) using 
a Chromel-Constantan thermocouple in a '*He cryostat 
for taking zero-field data above 4 K It is worthy 

of emphasizing that a careful pre-calibration of resistor 
sensors is indispensable for the precise thermal conduc- 
tivity measurements in high magnetic fields and at low 
temperatures. 

Figure 2 shows the temperature dependences of and 
Kab in zero field and several magnetic fields up to 14 
T. Like in usual insulating crystals there is a clear 
phonon peak at 8-9 K in both Kc and Kab, for which the 
peak magnitude is 8.5 and 24 W/Km, respectively. It 
can be seen that the heat conductivity of this organic 
crystal is rather large compared to common organic ma- 
terials; actually, the phonon peak in DTN is compara- 
ble or even larger than that in many inorganic crystals, 
like transition-metal oxides. One remarkable behavior of 
k{T) in zero field is a "shoulder" -like feature at ~ 2K; 
in addition, the "shoulder" moves to lower temperature 
upon increasing the magnetic field. This kind of temper- 
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FIG. 2: (color online) Temperature dependences of thermal 
conductivities of DTN single crystals for and Kah in both 
the zero field and several different magnetic fields up to 14 T, 
which applied along the c axis or the ah plane. The dashed 
lines in panels (b) and (d) indicate that the zero-field thermal 
conductivities at subKelvin temperatures show a T^ '' tem- 
perature dependence. 



ature dependence in niT) usually indicates a resonant 
phonon scattering f^Hl by some lattice defects, magnetic 
impurities, or magnon excitations, etc. Apparently, the 
sensitivity of the "shoulder" to the applied field suggests 
the magnetic origin of this resonant scattering. Another 
important feature of k(T) is that at subKelvin temper- 
ature regime, the thermal conductivity show a T^'^ de- 
pendence, which is close to the law of the standard 
temperature dependence of the phonon thermal conduc- 
tivity in the boundary scattering limit [25| . One possible 
reason of the slight deviations from law is due to the 
phonon specular reflections at the sample surface 26, 2^. 

Detailed magnetic-field dependence of the low- 
temperature thermal conductivity is a key to under- 
standing the mechanism of heat transport in the low- 
dimensional spin systems [2^. Figure 3 shows the k,{H^ 
isotherms for both Kc and Kab\ in each case the mag- 
netic field is applied both along the c axis and along the 
ah plane. Although the k,(H) behaviors in general are 
rather complicated in these four measurement configura- 
tions, one can easily notice the most striking result shown 
in Fig. 3(a), that is, at very low temperatures Kc{H) dis- 
play two peak-like anomalies across ~ 2.5 T and 12 T 
(11 c), which are very close to the reported critical fields 
Hd and Hc2 [13, [3, [I3| , and the anomalies are getting 
enhanced upon temperature approaching zero. Further- 
more, Kab{H) also show steep changes across these two 
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FIG. 3: (color online) Magnetic-field dependences of thermal 
conductivities of DTN single crystals at low temperatures. 



characteristic fields for H \\ c. The interesting point is 
that, both Kc{H) and Kab{H) do not show any drastic 
change at - 2.5 T and 12 T for 7? _L c. Since DTN 
does not exhibit any field-induced magnetic ordering or 
magnon BEC state when the magnetic field is perpen- 
dicular to the c axis 17[ , these data clearly indicate that 
the strong peak-like anomalies in Fig. 3(a) are related to 
the quantum phase transitions at Hd and i?c2- 

It is interesting to compare Kc{H) and Kab{H) behav- 
iors in the c-axis field, in which the magnon BEC can 
occur. Above 1.4 K, Kc and Hab have essentially similar 
magnetic-field dependences, as shown in Figs. 3(a) and 
3(c), while the difference between them shows up and 
becomes larger upon lowering temperature. For clarity, 
Fig. 4 shows a direct comparison of Kc{H) and Hab{H) 
at subKelvin temperatures. At 0.97 K, both Kc{H) and 
Kab{H) show a "U"-shaped curve: a steep decrease across 
2.5 T, a strong suppression but weak field dependence 
in the intermediate field regime, and a steep recovery 
of conductivity across 12 T. There is also a small dif- 
ference between these two curves, that is, a small and 
broad peak below 2 T shows up in Kc{H) isotherm. With 
lowering temperature, the behavior of Kab{H) does not 
change much, except that the suppression of the thermal 
conductivity in the intermediate field regime is gradually 
getting weaker and two shallow "dips" appear at 3 
T and 11.5 T. In the meantime, the behavior of Kc{H) 
changes much more drastically. First, the large peak- 
like anomalies at ~ 2.5 T and 12 T show up below 0.7 
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FIG. 4: (color online) Comparison of Kc{H) and K.ab{H) 
isotherms for H \\ c at subKelvin temperatures, (a) T — 
0.97 K. (b) T = 0.7 K. (c) T = 0.55 K. (d) T = 0.38 K. 



K and become more significant upon T ^ K. Second, 
the suppression of thermal conductivity in the interme- 
diate field regime is getting weaker rather rapidly with 
lowering temperature and finally evolves to an enhance- 
ment (compared to the zero-field conductivity) at 0.38 
K, which strongly suggests that there are two competing 
impacts on thermal conductivity induced by the mag- 
netic field. Apparently, the main difference between Kc 
and Kab can only come from anisotropic magnetic contri- 
butions to heat transport, acting as either heat carriers 
or phonon scatterers. Because of the strong anisotropy 
of the magnon dispersion 17|, it is a natural conclusion 
that the strong suppression of Kab{H) at Hd < H < Hc2 
is mainly due to the phonon scattering by magnons; on 
the other hand, although the magnon scattering can also 
weaken the phonon thermal conductivity along the c axis, 
the magnons (with stronger dispersion in this direction) 
can act as heat carriers and make an additional contribu- 
tion to the heat transport. Furthermore, between the two 
competing roles of magnons in affecting Hc, i-C, scatter- 
ing phonons or carrying heat, the latter one is apparently 
dominant at low temperatures. Note that because of the 
two competing effects of magnons on the heat transport. 
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the ability of magnons to carry heat must be much larger 
than what the raw Kc{H) data [H \\ c) demonstrate. To 
our knowledge, there has been no such clear evidence 
showing a large thermal conductivity in the magnon BEC 
state. 

Besides the above clear information demonstrated by 
the anisotropic heat transport behaviors, one may no- 
tice that the details of the k{H) data are actually rather 
complicated and some considerable further investigations 
are needed. In principal, the competing roles of magnons 
acting as phonon scatterers and as heat carriers may lead 
to complicated field dependence of k [H \\ c), including 
the peak- like anomalies of Kc{H) at Hd and Hc2 and the 
local maximum between two peaks. Besides, the peak- 
like anomalies can be closely related to the maximized 
magnon population at Hd and Hc2 , since the dispersion 
becomes quadratic at the critical fields while it is linear 
for Hci < H < Hc2 [i^ ; another contribution to the non- 
monotonic field dependence of Kc{H) between Hd and 
Hc2 may be coming from some upper magnon branches 
having rather small gap at the intermediate field 2^ . On 
the other hand, the Kc{H) and Kab{H) for _L c are also 
rather complicated. In general, the field dependences of k 
at very low temperature are rather weak, consistent with 
the fact that for _L c there is no field-induced transition 



and the magnetic excitations are always gapped [17|, |20| . 
The most drastic field dependence is the strong increase 
of Kab and Kc at 4-5 K, which is actually very similar to 
the behaviors in |j c. This is probably because in such 
high temperature region, where the temperature scale is 
comparable to the spin anisotropy gap, there are strong 
phonon scattering caused by magnetic excitations (re- 
member the resonant scattering feature of k{T) in zero 
field shown in Fig. 2), which can be weakened when ap- 
plied field increases the energy of magnetic excitations 

It is intriguing to point out that the experimental phe- 
nomenon cannot be simply explained as the strong mag- 
netic heat transport found in some low-dimensional spin 
systems [s^. Actually, it was originally predicted that 
the spin transport is diffusive and finite in the spin-1 
chain material since it is not an integrable system 3ll |. 
while the experimental results were quite controversial, 
both low and high magnon thermal conductivities have 
been observed in different compounds [s" 
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It 

is notable that in the magnon BEC state, the spin sys- 
tem is no longer one dimensional; instead, it is a three- 
dimensional ordered state [13, EM Because of the 
crossover of the dimensionality and the character of mag- 
netic quasiparticles at the BEC transition (ssj ] , it is natu- 
ral to expect a different mechanism of magnon heat trans- 
port from that of the low-dimensional systems when the 
field-induced long-range magnetic order is established. 

To summarize the main finding of this work, the 
magnon heat transport is found to be large in the magnon 
BEC state of DTN single crystals, pointing to a direct 



analogy between the magnon BEC and the conventional 
one. Since the BEC condensate does not carry entropy, 
the large low-T magnon heat transport can only be re- 
lated to the uncondensed part of magnons. An elabo- 
rate theory, probably based on the two-fluid model es- 
tablished for superfluid "^He [l5[, is called for quantita- 
tively describing the heat transport in this novel state of 
quantum magnet. 
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